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ABSTRACT 

Context. Galactic synchrotron emission, rotation measure (RM) and the deflection of ultra-high-energy-cosmic-rays (UHECR) permit 
detailed studies of the Galactic magnetic field (GMF). The synchrotron emission has also to be measured in order to be separated 
from other astrophysicaly interesting signals like the CMB. 

Aims. We present a publicly available code called hammurabi for generating mock polarized observations of Galactic synchrotron 
emission for telescopes like LOFAR, SKA, Planck and WMAP, based on model inputs for the GMF, the cosmic-ray density distribu- 
tion and the thermal electron density. We also present mock UHECR deflection measure (UDM) maps based on model inputs for the 
GMF. In future, when UHECR sources are identified, this will allow us to define UDM as a GMF probe in a similar way as polarized 
radio sources permit us to define rotation measures. 

Methods. To demonstrate the code's abilities mock observations are compared to real data as a means to constrain the input parame- 
ters of our simulations with a focus on large-scale magnetic field properties. 

Results. The Galactic magnetic field models in the literature seem to fail to reproduce any additional observational data which was 
not included in their design. 

Conclusions. As expected, attempts at trying to model the synchrotron, UHECR deflection and RM input parameters, show that any 
additional observational data set greatly increases the constraints on the models. The hammurabi code addresses this by allowing to 
perform simulations of several different data sets simultaneously, providing the means for a more reliable constraint of the magnetized 
inter-stellar-medium. 

Key words. Radio continuum: ISM ~ ISM: mangeUc fields - ISM: cosmic rays 

1. Introduction ing a publicly available softwarq] capable of generating mock 

synchrotron and Faraday rotation observations as well as mock 

There are several different observations which probe the Galactic UHECR deflection maps. 

magnetic field (GMF): synchrotron radiation, rotation measure, 

c3 UHECR deflection, dust related observations (e.g. starlight po- The code is constructed in such a fashion that it should be 

■ ■ larization iHdlii (i2000) . polarized dust emission), and more feasible to adapt it to do any sort of all or partial-sky line-of- 

recently a tom alignment spectro scopic observations were pro- sight integral observations, (see, e.g., [Waelkens et al.||2008| and 

posed (see l^n & Lazarianf2006l) . | Sun et al. |i 2008 , ) . 

A better understanding of the inter-stellar magnetized We wiU give first a brief account of previous work on GMF 

plasma and the radio emission processes related to it is also modeling in Sect. |2] and then we will give a description of the 

paramou nt for current and future CMB experim ents (see, e.g., involved Faraday rotation, UHECR deflection and synchrotron 

[Mivifle- Deschenes et al.ll2008t iPage et al1l2007l) . where in par- physics, and the simplifying assumptions used in our computa- 

ticular measurements of the polarized CMB signal wiU be lim- tion (Sect. |3). The implementation of a line-of-sight (hereafter 

ited b y our knowledge of foreground emission dTucci et all LOS) integration scheme which mimics the radiative transfer is 

l2005h . discussed in Sect. |4] We present some test output using simple 

Observational studies of the GMF are usually based on one standard magnetic field, thermal- and cosmic-ray electron mod- 
single sort of observable. However, since the diff^erent observ- els. We disc uss the role of the turbulent field and sub-grid model- 
ables and measurement methods provide complementary infor- i"g i" Sect. |5.1.4| and the one of helical magnetic fields in Sect, 
mation, constraining the GMF is logicaUy better if aU possi- EES] Finally we conclude in Sect.|6] 
ble observations are considered simultaneously. We address the 
need to confront those observations with models by present- 
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Software available at 



email: waelkens@mpa-garching.mpg.de |http : //www . mpa-garching . mpg . de/~waelkens/haitimurabi | 
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2. Previous work on Galactic magnetic field 
modeling 

GMF modeling already counts with several c ontrib u tions from 
the communi ty. See , e.g., Beu ermann et al. l (Il985h: iHan et al 



(2006); Bro wn et all (120 03); Ha verkorn et a l. (2006); lPage et al. 



62007i) ; iSun et al.l (l2008i) and many others. The data used are 
mostly a number of total and polarized emission radio sur- 
veys of the Galaxy dJaylor et al. 2003; Gaensler et al. 200J| 
iHaslam et al"] 119821: iReich & Reich 1986.; Page et al.. .2007.) as 
well a s measurements of Faraday rotation (iDineeri & ColesI 



120051; iJohnston-Holhtt et all 120041: iHan et al.1 l2006l) and 



starlight polarization catalogue compiled by iHeilesI (l2000h . 
Theoretical predictions are also heavily based on experience ob- 
tained f rom synchrotron observations of other spiral Galaxies 
(see e.g. lBeck et al.lll996h . Our work complements previous ef- 
forts by making it possible to compare simultaneously the largest 
possible number of observables to theoretical predictions. 



order of a few micro Gauss with a yet unknown topology. In the 
following, a couple of simplifying assumptions are made: 

- The relativistic CR electrons have an isotropic velocity dis- 
tribution, as is measured with high accuracy to be the case 
at our location in the Galaxy (their propagation, though, is 
not isotropic due to th e presence of ordered magnetic fields; 
lYan & Lazarianll2008l) . 

- The cosmic-ray electron spectrum is assumed to be a power 
law with spectral index p. This widely used simplifica- 
tion is motivated by the theory of shock acceleration (a.k.a. 
Fermi a cceleration, which predicts power law energy distri- 
butions; iDrurvl 11983.) and is simultaneously confirmed by 
the measured cosmic-r ay electron spectrum at Earth (e.g. 
iGaisser & StanevI l2004 and references therein). The same 
observations, as well as sophisticated simulations have, how- 
ever, shown that the power law as sumption is not ad equate 
for the entire energy spectrum (see lStrong et al.ll2007h . 



3. Physics included in the code 

Here we describe the physical processes included and the as- 
sumptions underlying the code. 



3.1. Faraday rotation 

The polarization angle of an electromagnetic wave is rotated 
when crossing a magne tized plasma. This effect is known as 
Faraday rotation (see e.g. lRvbicki & Lightmanid 19791) ). The ob- 
served polarization angle ;^' will be a function of the rotation un- 
dergone when crossing the magnetized plasma, the square of the 
observation wavelength and its original (or intrinsic) angle xo at 
the polarized source. 



X^RMA^+Xi). 



(1) 



The rotation measure (hereafter RM), which quantifies the hnear 
rate of change of the angle ;t' as a function of A^, is a function of 
the integral of the magnetic field B^os along the LOS weighted 
by the thermal electron density n,,, 

pthere 

RM = ao dr He Bios, (2) 



where oq — q^jilmny ), irie is the electron mass, c is the speed 
of light, and qe is the electron charge. 

The RM can be measured directly via a fit to Eq. [1] only 
in the particular case of a Faraday screen, i.e., if the observer 
and the polarized source have a cloud of magnetized plasma in 
between them, but no source-intrinsic Faraday rotation occurs. 
In the more complicated scenario of several sources along the 
LOS embedded in the magnetized plasma, the RM cannot be 
measured in that way since the polarization angle will no longer 
obey a linear dependence onA^.A complicated frequency depen- 
dence of the polarization angles arises in such a case, and that 
is the situation which is typically found in our Galaxy. There, 
the synchrotron-emitting cosmic-ray electron population is em- 
bedded in the magnetized plasma that is producing the Faraday 
rotation simultaneously. 



3.2. Synchrotron emission, total and polarized 

Relativistically moving charges in a magnetic field emit syn- 
chrotron radiation. In the Galactic case we deal with a cosmic- 
ray electron population mostly arising from supernova explo- 
sions and subsequent shock acceleration, and the GMF of the 



The synchrotron emissivity (i.e. power per unit volume per 
frequency per solid angle) is partially linearly polarized. Its in- 
tensity and polarization properties depend on the strength and 
orientation of the perpendicular (with respect to the LOS) com- 
ponent of the magnetic field, Bj^, and the cosmic-ray electron 
spatial and energetic distribution. The emissivities are usually 
subdivided into two components, 7j_j| - dEj_\\/dtda>dQ.dV, 
respectively perpendicular and parallel to B^ ps , following 
iRvbicki & LightmanI (IT979I) and lWestfoldl (fT959l) : 
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Here C depends on the position in the Galaxy and is defined by 
N{j)dy - Cy^''dj, y being the Lorenz factor, N(y) the number 
density of electrons e [y,y + dy], and p is the spectral index as 
mentioned above. The charge of the electron is given by qe, the 
mass by m^. and w - 2nv, where v is the observation frequency. 
The specific intensity / as a function of observation frequency 
and LOS direction h is 



I{io, h) - \ dr \ji_(cx), rh) +j\\(cx), rh\ , 
Jo 



(5) 



and the polari zed specific intensity P expressed as a complex 
variable is (see lBurnlll966l) : 



-lixirii) 



dr (j^(co,rh) -j\\{oj,rh))e 





The intrinsic emission polarization angle X() is given by 
the inclination towards the Galactic North Pole of the local 
perpendicular-to-the-LOS component of the magne tic field at 
each position in space (same convention as adopted in lPage et al.l 
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l2007h . The Stokes I, Q, U parameters^ are then the integrals over 
sohd angle Q.: 



I' 



I = I dm, 

and 

Q + /U= I dQP. 

3.3. UHECR propagation 

UHECR's are deflected by the GMF due to the Lorentz force. 
The Larmour radius for relativistic particles is rg - p c/Z qg Bj., 
Z qe being the charge of the UHECR, and p the momentum per- 
pendicular to the magnetic field Bj_. Provided, the sources of 
UHECRs could be identified, an UHECR deflection measure 
(UDM) can be extracted from the UHECR arrival distribution 
via fitting the arrival data (position & energy). The net deflec- 
tion of an UHECR can be app roximated by the LOS integral 
(see e.g. lKachelrieB et a~ 



I be app i 
III2OO7I) 



©offset ^ ffl'/r-i = ^ {dlB^ = ^UDM. 
J ' pc J pc 



(6) 



4. Implementation 

In this section we present the implementation of the physics de- 
scribed in Sect. l3.2l and the technical characteristics of the ham- 
murabi code. Given a 

- 3D GMF model, 

- 3D cosmic -ray electron density model, 

- 3D thermal electron density model, 

HAMMURABI computcs full sky maps for 

- the Galactic RM contribution to the extra-galactic sources, 

- synchrotron I,Q and U Stokes parameters, taking into ac- 
count intrinsic Galactic Faraday depolarization and, 

- UHECR deflection intensity and orientation maps. 

The sky maps in hammurabi are subdivided into equal 
area pixels f ollowi ng the HEALPi?o pixelization scheme of 
iGorski et al.l (l2005b . The total number of pixels for an all-sky 
map, which defines the angular resolution, is A^pix = 12NSIDE^, 
with NSIDE = 2*^ and ;t = 0, . . . , 13 (a HEALPix package Hmi- 
tation, which can be altered, see App.lAt. The angular size of a 
pixel (A0) can be approximated as 



A9: 



3 3600' 
n NSIDE 



(7) 



The observation volume associated with one of the 
HEALPix-map pixels is sampled along the LOS at a constant ra- 
dial interval Ar. As a consequence of the conical shape of the im- 
plied effective observation beam, the volume units increase with 
radius, and hence the weights of the sampling points along the 
LOS do also increase. To limit the amount of non-homogeneity 
of the sampling, the code allows the volume resolution to be 
increased by splitting the beam in sub-beams at some radius. 



^ Here the specific intensities are in itahc (I,Q,U), while the observed 
Stokes parameters (I,Q,U) are denoted by Roman letters. 
^ http://healpix.jpl.nasa.gov 



which subsequently can be further split later on. We call the im- 
plied 3-D sampling grid the "3D HEALPix grid". See further 
details in App.lAl 

Formally, the maximally achievable volume resolution, i.e. 
the largest volume unit at the finest possible resolution A0 ~ 
0.43', (corresponding to NSIDE = 2'^) is y„, = 47rRl„^Ar(12 ■ 
2^*)"' ~ (4pc)^ for R,„a_^ ~ 32kpc and Ar - 4pc since we im- 
pose an approximately cubic volume unit according to Ar ~ 
(4;r/?^„^/12 • 2^*)'^^. This implies that variations of the input 
parameters (magnetic field, thermal electron distribution,...) on 
scales smaller than this volume can only be taken into account 
with sub-grid modeling if they persist to the largest radii simu- 
lated. Note however that this is an upper limit, given that fluctu- 
ations are likely to be stronger closer to the Galactic disk, where 
resolution, due to the cone-like shape of the observation volume, 
is going to be higher anyways. 

The i ntegration is per formed assuming an optically thin 
medium. I Sun et al] (l2008l) enhanced the code by introducing 
free-free absorption. This implementation, relevant mainly at 
frequencies below a GHz, is present in the code but not described 
here. Furthermore. Sun et al.. (.2008.) also introduced a coupling 
between the thermal electron density and the random component 
of the RM. This implementation is necessary to explain the de- 
gree of depolarization at thek simulated 1 .4GHz map. It is also 
present in the c ode but not described here. For further details see 
ISun et all ( I2OO8 ). 

The 3D HEALPix grid, since the value of one observation 
pixel stems from the contribution of several sub-beams, allows 
one approximately to take into account effects like beam depo- 
larization. See App. lA.ll for details. 

4.1. Features 

HAMMURABI is also suited for simulations of partial sky cover- 
age. A single pixel or a list of pixels representing either patches 
or separate locations on the sky output maps can be computed. 
An extension to that is the option to compute RMs for an indi- 
vidual pixel only up to a certain specified distance. This is rel- 
evant for simulating RM observations of radio pulsars or other 
polarized sources with distance information in our own Galaxy. 
Polarized sources along the same LOS allow us to probe frac- 
tions of the ISM plasma, unlike extra-galactic polarized sources 
which give us the integrated RM along the entire LOS through 
the Galaxy. The mock observations for this case are done without 
beam-splitting and at highest possible angular resolution, since 
RM are effectively obtained from point sources. 



5. Examples 

In this section we choose standard 3D input models for a demon- 
stration of the code's results. It is not the scope of this work to 
present new scientific findings, but to present a proof of concept 
of what can be done with hammurabi. An actual application of 
the code has been done in lSun et all (l2008h . and we will briefly 
refer to some of their results here. Their preferred GMF model 
and cosmic-ray electron model, derived by trying to fit a broad 
range of observ ables, are c ompared to the corresponding mod- 
els suggested bv lPage et al.l (12007.) . fitted to reproduce solely the 
polarization angles observed by the W MAP sat ellite. Note tha t 
it has been shown bv lSun et all ( l2008h that the Page et all (l2007h 
model, as well as every model from the literature analyzed in 
their work, fail to reproduce different sorts of observations of 
the magnetized ISM satisfactorily, since they are all constructed 
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by considering only a particular sort of data. The discrepancies 
of the models presented here, visible by eye, are suitable for dis- 
playing the codes abilities. 



The same simulations also present a seizably more sophisti- 
cated cosmic-ray electron distribution. Incorporating the output 
of those simulations into hammurabi is planned as an extension. 



5.1. Inputs 

5.1 .1 . The thermal electron density model 

For t he therm al electron distribution we use the NE2001 model 
(see ICordes & Lazio 2002) . This model subdivides the Galaxy 
into several large-scale structure elements like a thin disk, a thick 
disk, spiral arms, as well as some local small-scale elements such 
as supernova bubbles. 



5.1.2. The cosmic-ray electron density model 

We use the cosmic -ray electron density models suggested by 
iPage et al. (2007) and Sun et al. (2008). 



- The model in iPage et all (l2007h consists of an exponentially 
decaying disc with characteristic hight h^ = 1 kpc and a char- 
acteristic radius of hr = 5 kpc. Note that to compute the 
synchrotron emissivity, as can be seen in Eqs. |3]and|4] we 
need the spatially dependent function C, not the cosmic -ray 
electron density. However, since we are assuming a Galaxy- 
wide unique power-law energy spectral slope with a spec- 
tral index p = 3 (note that in principle the code allows to 
associate a different spectral index for each volume unit), 
these quantities are proportional to each other. The value 
of CEaith - 6.4 ■ lO'^cm'-^ at Earth 's position is observed 
(see e.g. Fig. 4 of IStrong et al.ll2007l) . and although it is not 
clear that it is representative for other regions in the Galaxy 
dSun et al . 2008; Pohl & Esposito 1998; Strong et al. 2004), 
we use it as a zeroth order approximation f or the normaliza - 
tion (which is not necessary/provided in .Page et al.l (l2007h ) 
of our distribution. 



C - Co exp [-r/hr] sech (z/hd) . 



(8) 



Here r is the Galactic radius, while z is the hight. Co is such 
that C - C Earth at Earths position. 
- ISun et alj (i2008) proposes 



C(R,z)^CQexp 



R-Rf^ 



8 kpc 1 kpc 



(9) 



with Co = CEarth, whilc C{R < 3 kpc) = C(R = 3 kpc) and 
C(\z\ > 1) = 0. The abrupt truncation at a scale hight of 
|z| > 1 as discussed by Sun et al. (2008) is necessary to aco- 
modate a low synchrotron emission at high latitudes where 
an anomalously strong halo magneti c field is re quired to ac- 
count for high RM measurements. I Sun et al] (2008) warn 
that this seems unrealistic and suggests that a larger scale 
hight of the thermal electron density could resolve the prob- 
lem by diminishing the required value of the halo magnetic 
field. They conclude that the high-latitude thermal electron 
density should be better investigated. Furthermore the spec- 
tral index p is 3 for obs ervation freque ncies larger than 408 
MHz and 2 otherwise. ISun et al.l ( I2008.) adopt t his si mpli- 



fication based on obse rvations of iRoger et al.l d 19991) and 
iReich & ReichI (Il988allbl) . which support a flatter spectrum 
below 408 MHz. 

As mentioned in Sect. 13.21 simulations (see the GALPROP 
code bv lStrong et alj2007l) show that the assumption of a power- 
law energy distribution is not applicable to the entire spectrum. 



5.1.3. Large-scale Galactic magnetic field models 

The behavior of the synchrotron emissivity is mainly driven by 
the magnetic field distribution, as shown in Eqs. [3] and H] 

It is a common practice to subdivide the Galactic magnetic 
field in a large-scale and a small-scale component. The latter will 
be addressed in the next section [5.1. 41 The subdivision in these 
two classes of fields is somewhat arbitrary, and we adopt here the 
convention that the large-scale field is statistically anisotropic at 
any point in the galaxy, while the small-scale field is not. 

GMF modeling has b een done in direct st udies of the galac - 
tic magnetic field. Se e IFage et alj (l2007l); ISun et al.1 (l2008l); 
Jansson eFal] (|2007|); [Haverkorn et all (l2006l) ; iBrown et al.1 
(i2003i) ; iHan et al.1 (12006'), and othe rs. In the context of UHECR 
propagation (see IJCachelrieB et al] 120071; iTinvakov & TkachevI 
llOOllProuza & SmidalMBblHarari et al.l2000h . Typically GMF 
models, inspired by observations of the radio-polarization pat- 
terns of other spiral galaxies, present spiral-U ke structures. Here 
we present simulations using the model from Page et al. (2007h 
and th e favored (ASS+RING) model proposed by Sun et al] 
(|2008|) . The expressions are reproduced in the App. |B] Further 
magnetic field models can easily be incorporated in the code. 



5.1.4. Small-scale magnetic field 

Although the turbulent component of the GMF has been stud- 
ied extensively (see e.g. 'Haverko rn et al1l2003l ; iHan et al]|2004l ; 
Haverkorn et al. 2006, 2008) there is, to our knowledge, not 
yet any model which describes it to a satisfactory degree. 
Observational constraints on its relative strength compared to 
that of the large-scale field put it at roughly t he same ampli- 
tude (i.e. a couple of micro Gauss; lBeckll200lT) . hammurabi al- 
lows one to simulate, given a magnetic field power spectrum, 
a Gaussian random field realization. A realistic turbulent field, 
however, needs not be Gaussian, since simulations of turbulence 
show sheet like structures evidently different from what is ex- 
pected from a Gaussian field (see e.g. lSchekochihin et al.ll2004l ; 
ISchekochihin & Cowlev..2006.) . This means that the Gaussian 
random field should per construction reproduce the correct two 
point statistics of a real field, however it will probably not repro- 
duce higher order statistics as r eal fields are not observed to be 
Gaussian (see e.g. the figures in lJaegerslll987l ; IClarke & EnsslinI 
2006) . Furthermore, this modeling assumes the power spectrum 
of the field to be known, which in the case of the Galaxy, to 
our knowledge, has not yet been determined. There are, how- 
ever, measurements o f the magnetic power spectrum in the intra- 
cluster medium (e g. IVogt & EnBlinll2005l: iGuidetti et ai]|2008l; 
iGovoni et al]|2006l) . 



5.1.5. Helical fields 

Helical magnetic fields have a broken mirror symmetry and can 
best be visualized by twisted flux tubes. Topological helicity is 
given by J dV A ■ B, where A is the vector potential of the field, 

B - V X A. Current helicity is given by j dV j ■ B, where j is the 
electrical current, satisfying j -VxB. Both types of helicity are 
closely related. Topological helicity is a nearly conserved quan- 
tity, even in non-ideal MHD, and plays a crucial role in mag- 
netic dynamo theory. Therefore hammurabi was adapted to read 
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in random small-scale fields with magnetic helicity generated by 
a separate code. This code, also available for download, is called 
GARFiELDS (first applied in Kitaura & Ensslin 2007), where helic- 
ity is imprinted onto a divergence free Gaussian random field B' 
by 

Bm = -^=^= {6jm + iaeijjci^ B'j. (10) 



polarised intensity, non-helical field 



vr 



+ a^ 



Here B'- = Bl(k) and it is such that (B'-Bp = {dij - kikj) P(k)/2, 
where P{k) is the magnetic power spectrum. The hat denotes a 
Fourier transform, 6 is the Kronecker delta, k, is an unit vector 
and eijk is the permutation symbol. Maximal helicity is found for 
a = -I and 1. This leads to a magnetic correlation tensor of the 
form 



(Mj) 



-[6u-kikj)-ieij,k,-^^ 



P(k). 



(11) 



5.2. Synchrotron and RM outputs: comparison with 
observations 

In this section we compare our example simulation with ob- 
served data. The mock total intensity I, polarized intensity P, 
polarization angle PA and RM maps for our two diff'erent field 
configuration s and s lightly different cosmi c-ray e lectron distri- 
butions of Sun et al.l (12008) and Page et al] (l2007h are shown in 
comparison to the observed synchrotron emission at 0.408 MHz 
(total intensity) and 23 GHz (polarized emission), and a RM map 
extrapolated from the existing still sparse RM observations (Fig. 
[1). All maps in this figure and throughout the paper are at a res- 
olution of NSIDE = 128, Ar = 0.21kpc, i/obs = 2, i/,™x = 3 
(see definitions in App. |A|. A Gaussian random field with a 
Kolmogorov spectrum k^^^^ with a lower cutoff at ko - Ikpc"' 
(see, e.g., Han et al. 2004), and an upper cutoff given by the 
Nyquist frequency of the simulation box kffyquisr = 3.2kpc"', 
was added on both large-scale GMF models. We caution that 
this scale range is insuflicient to describe the turbulent field in 
the Galaxy, as can be seen when noticing the even smaller-scale 
fluctuations present only in the observ ations, however it serves 
our demonstration purposes. Following lSun et al](l2008l) the ran- 
dom field amplitude is (B^„„)'''^ = 3;uG. 

Magnetic field models in the literatur e are mostly designed 
to fit a selected set of data. For example iPage et all (l2007h fit 
PA observations, and, as can be seen in Fig. [1] does this re- 
markably well. However, their model evident ly fails to repro - 
duce any of t he other ob s ervati ons (as shown bv lSun et al.l2008l) . 
Furthermore I Sun et al] (l2008h note that most spiral magnetic 
field configurations have rather similar PA, suggesting a highly 
degener ate fitting problem . 

The ISun et alj (|2008|) models for the GMF and the cosmic- 
ray electron distribution were designed to fit longitude and 
latitude RM profile observations and I and P observations at 
408 MHz, 1.4 GHz and 22.8 GHz, but no optimization was per - 
forme d for the PA observations. Note also that . Sun et alj (|2008|) . 
unlike iPage et all ( l2007l) . who did a parameter scan, adjust the 
model parameters based on trial and error. The attentive reader 
might notice a faint horizontal stripe like pattern seen in the I 
and P simulations of this model, which is due to a combination 
of our limited radial step size Ar = 0.21kpc, and the abrupt cut of 
the cosmic-ray electron density at 1 kpc away from the Galactic 
plane (as discussed in Sec. |5.1.2] i. In other words, the code can- 
not resolve a sharp cut. A smaller Ar, of course, makes the fea- 
ture vanish, though it would be inconsistent with our recommen- 
dation in App.|A]to peg Ar to the NSIDE parameter. Probably the 




polarised intensity; helical field 




Fig. 2. Polarized intensity of a field without and one with he- 
licity. The colors in the plot are equally distributed over area to 
enhance the visualization of structures at the expense of propor- 
tionality information. 



most striking diff'erence between the two models lies in the dif- 
ferent RM maps , mainly due to so called field reversals in the 
ISun et al.l (l2008h model, which are absent in the WMAP model. 



5.3. Helicity in the ISM? 

Our code provides the possibility to study the effects of mag- 
netic helicity on observations. In Fig.|2]a field realization with 
and without helicity are compared. Although they are clearly 
different, there is no qualitative difference between the single 
frequency maps from fields with helicity and without. These re- 
sults reinforce that measuring helicity in the ISM is difficult and 
the necessary multi frequency meth ods to extract it e.g. from 
Faraday tomography observations (B rentiens & de B ruvn 20051: 
:Schnitzeler et al...2007.) still have to be developed. The availabil- 
ity of a tool for generating mock observations might be of as- 
sistance for the development of such methods and for feasibility 
studies of observations aiming at determining the helicity in the 
ISM. 



5.4. UHECR deflection outputs 

HAMMURABI as a generic tool to study the GMF also helps to pre- 
pare for UHECR-based magnetic field studies, once the sources 
of the UHECR particles are identified. We included the option to 
compute the deflection of ultra-relativistic charges by the GMF 
in the code. This is yet another observable by which the GMF 
may be further constrained in the future. Presently there are only 
a couple of hundred registered events and their origin is still 
speculative (iThe Pi erre Aug er Collaboration et alj|2007h. Fig. [3l 
repres ents the deflection intensity for the large-scale ISun et al] 
( |2008|) model. Note that the strong halo field is controversial (as 
discussed in Sect. 15.1.21 ). The deflection angle of an individual 
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Fig. 1. Comparison of mock observations based on models from lPage et al.1 (l2007h and lSun et alJ (l2008h with observed data. Each 
row is respectively total in tensity I, polarized intensity P, polarization angle PA and rotation measure RM. The third column contain s 

(1 198 111 1982) . P and PA are WMAP observatio ns (iPaee et al.ll2007t iHinshaw et alJl2008l) , 



the observables. I is from Haslam et al. 



2005h using the lSimard-Normandin et all (Il98lh data. 



UHECR is obtained by multiplying the deflection intensity map 
with Zqe/pc of this particle, as in Eq. |6] 



5.4.1 . A LOFAR, SKA, WMAP, Planck tool 

The frequency range of hammurabi synchrotron simulations, 
from ~ 100 MHz up to ~ 100 GHz, covers currently running 
experiments like WMAP as well as the upcoming generation of 
radio telescopes: Planck, LOFAR, SKA and the like, hammurabi 
has been developed to support the scientific exploitation of the 
data of these experiments and to provide synthetic observations 
for design studies. It is capable of generating full- or partial- 
sky (not shown) maps as well as individual LOS (see Figure [S]), 
which are useful for non diffuse measurements, e.g. RMs. 

Figure |4] shows the full sky polarized intensity P for the 
aforementioned Sun et al. (2008) model (applying their low fre- 
quency corrections and coupling between thermal electrons and 
the random RM component). 




Fig. 3. UHECR deflection measure map from the I Sun et al] 
( |2008|) large-scale magnetic field model. Currently only a couple 
of hundred UHECR detections have been reported. 
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Fig. 4. Simulated polarized intensity P at 1 and at 2GHz of the 
I Sun et al] (l2008l) model. The colors in the plot are equally dis- 
tributed over area to enhance the visualization of structures at the 
expense of proportionality information. The Faraday depolariza- 
tion effect changes the complex patterns as a function of the 
wavelength. The higher thermal electron density on the Galactic 
plane keeps Faraday depolarization active long after it already 
ceased to be relevant in the Galactic halo. 




Fig. 5. The simulated rotation measure as a function of disper- 
sion measure. Here we used a Gaussian random field with an 
uppercutoff of 5 kpc and a root-mean-square of 6 //Gauss. Pulsar 
polarization measurements are, with the advent of SKA, ex- 
pected to become very numerous, allowing us to chart the GMF 
in three-dimensions, as is shown in the present plot. 



In principle, multi-frequency simulations permit one to take 
into account beam-width depolarization effects. Since the mea- 
surements of real telescopes correspond to an integral over 
some frequency band, band-width depolarization will also hap- 
pen provided the Faraday effects are strong enough across the 
band-width. Although band-width depolarization is not a stan- 
dard HAMMURABI feature, it can be emulated by coadding closely 
spaced frequency maps. 



The framework of the code is built t o allow addition s of an y 
other beam-like observables such as in IWaelkens et alj ( l2008l) . 
fSun et al., ( .2008,) . 



6. Conclusion 

We have presented the hammurabi code, a tool for simulat- 
ing rotation measure maps, total and polarized Galactic syn- 
chrotron emission maps taking into account Faraday rotational- 
depolarization effects, as well as UHECR deflection maps. 

Drawing a set of well known input models from the liter- 
ature (NE200 1 for the therm al electron distribution, I Sun et al] 
(I2OO81) and iPage et al.l ( |2007|) for the cosmic -ray electron distri- 
bution and the Galactic magnetic field), we show example out- 
puts comparing them to corresponding observations to illustrate 
the code's abilities as a scientific tool for charting the Galactic 
magnetized plasma and the cosmic -ray electron distribution. 

Full galaxy simulations are currently limited by the finite 
HEALPix-grid resolution, thus being unable to probe fluctu- 
ations with characteristic lengths of less than 4pc. A second 
caveat is that the widely used approximation of a power-law en- 
ergy distribution of the cosmic-ray electrons has been shown to 
be ques tionable by sophis ticated cosmic ray propagation simu- 
lations JStrong et al.ll2007l) . The degree to which this affects the 
precision of our simulations has yet to be assessed. Furthermore, 
our non-ray-tracing calculations of UHECR deflections might 
not fulfill strict precision requirements as soon as the UHECR's 
trajectory through the Galaxy is significantly bent; they serve 
rather as a first approximation. 

It was shown that, unsurprisingly, models designed to fit only 
a fraction of the available observational data on GMF might fail 
to reproduce the remaining observational information not taken 
into account in their construction. 

This is an indication that the constraints to the GMF might 
be highly degenerate. Hence confronting the models with the 
broadest possible range of observations on the Galactic magnetic 
field is paramount to achieving any useful statement about the 
field. 

Furthermore we have displayed Hammurabi's capabilities of 
generating low-frequency mock observations, where Faraday ef- 
fects play a significant role. This can be used for feasibility stud- 
ies and the analysis of the actual observations of forth coming 
low-frequency radio telescopes like LOFAR or the SKA (note 
that the latter can also observe at high-frequencies). 

HAMMURABI Can be applied for constraining the Galactic 
magnetized plasma and increasing our understanding on radio 
Galactic emission and UHECR source locations. Already a wor- 
thy scientific goal by itself, it will also have implications on our 
understanding on foreground emission subtraction for experi- 
ments like Planck or WMAP 
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Appendix A: The 3D HEALPix grid 

The HAMMURABI code generates HEALPix maps with an angu- 
lar resolution defined by the parameter NSIDE - 2*, where 
k e [0, 13] (the maximum k is a computational limit; it could 







obs 1 


______ 




n\ 


^^^^ 


^ 


observer 








W""^ 


Ar 







Fig. A.l. An example of the 3D HEALPix grid. For simplicity 
it is presented in 2D. The gray shaded area corresponds to the 
integration volume. 



be extended only by alterations to the original HEALPix pack- 
age). The integration volume has a cone-like geometry. To min- 
imize the non-homogeneous volume sampling induced by that, 
the integration volume is consecutively subdivided in so-called 
sub-beams. Each section of these sub-beams is contained in a 
shell centered on Earth. We shall call this grid the 3D HEALPix 
grid. 

This subdivision will depend on 

- the radial length Ar of each volume unit in the observation 
beam, 

- how many shells //max in total there will be, 

- what resolution, or equivalently what NSIDEobs the observa- 
tion shell has, 

- and finally the shell number of the observation shell //ob.s- 



See an example in Fig. lA.ll 

The distances at which the beam is split, or equivalently the 
upper boundariesQ of the shells, are given by 



J _ p Tf-Hmax 



(A.l) 



Where the index / e [l,//max]- Due to the disc-like shape of our 
galaxy, the eff'ective radius out to which we perform our compu- 
tations varies and decreases sharply towards high latitudes and 
towards the anti-center direction. 

The maximum volume for each shell is a constant corre- 
sponding approximately to 



V™ 



Ar 



= And], ^ T- 

""^ 12 ■ NSIDEi 



(A.2) 



We suggest a definition of Ar such that the maximum volume 
unit at the largest shell is close to a cubic form. 

Note that the maximum volume unit Vmax decreases with 
higher observation resolution. It defines the minimum scale 
which can be resolved with the code for a certain resolution. 
Fluctuations smaller than that cannot be resolved by the code. 
This is a formal upper limit, since the Galactic regions more 
prone to have small scale fluctuations are likely not at the far 
end of the observation volume, where V„jax is located, but instead 
at regions closer to the observer, where the spatial resolution is 
finer. 



^ The code will automatically round down di to a multiple of Ar. 
'' Choices of any other constant are of course possible. 
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(A.3) 



A. 1 . Line of sight integration 

For the i-th volume element of a cone, according to the formulae 
in the theoretical section |3] we compute the following, 

RMi - aoneiBi^ios^r 

Qi = PiCos{2xi) 

Ui = Pism(2xi) 

Bjc,i = Bicos(xi,o) 

By J = Bi&mixi^Q) 

Note that the calculated Stokes values {/, and Qi for each vol- 
ume element include the effect of foreground RM. The CR elec- 
trons are assumed to follow a power law with an energy spectral 
index p (see Sect. 13.21 ). The values C/ and Cp are dependent on 
the spatialy dependent spectral index p and the spatial distribu- 
tion component of the cosmic rays C, and can be obtained from 
the formulae given in section[3]and lRvbicki & Lightmanl ( ll979h . 
The intrinsic polarization angle ;t',o is locally defined as zero if 
the magnetic field is pointing in positive Galactic-coordinates l- 
direction and 7r/2 when the magnetic field is pointing in positive 
Galactic-coordinates /j-direction. 

The intensities, RMs and UDMs for a pixel are straightfor- 
ward to compute. They consist of an integral of the contributions 
from all the volume units as below (Eq. |A.4t . For the intensities, 
in case of a subdivision of the observation cone into sub-cones, 
the computation is done by averaging the set of sub-cones. This, 
however, is not done for the UDM and RM values, since they are 
respectively observed as being linearly dependent on A^ and A^, 
where A is the wavelength. The beam averaging could destroy 
that linear behaviorQ. The UHECR deflection measure is given 

by UDM = JuDMl + UDMj, while the deflection orientation 
is given by ©defi = arctan 



UDM,: 

UDMy ■ 



Q = Z/G, 

u = T.i Ui 

RM = 'LiRMi 

UDM, = Y,iB,,jAr 

UDMy = i:,B,,,Ar 



(A.4) 



- iPage et all (l2007h writes 

B(r,(f>,z) = Bo[cosi/'(r) cos;^(z)r-l- 

siniffir) cosx(z)^ + sin x(z)z\ 



(B.l) 
(B.2) 



Herei/'(r) = ipo + i//iln(r/8kpc), x(z) - Xo tanh(z/lkpc). The 
radial variable r € [3kpc, 8kp c1, yp = 25°, lAi = 0.9° and 
tf/Q - 27°. Bo is not specified in lPage et al.l (|2007|) and we put 

Bo = 4^G. 

While lSun et al.l(l2008l) presents some suggestions for fields, 
and we pick their ASS+RING parameterization: 



BP = Di{R,4>,z)D2iR,(l>,z)smp 
B'? = -Di(R,<f>,z)D2(R,cf>,z) cos p 
,B? = 
Where 



(B.3) 
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Di(r,z) 



Boexp 



Re, 



Ro 



Br 



-] r>Rc 

Zo) 

r<Rr 



(B.4) 



Here Rq - 10 kpc, zo = 1 kpc. Re = 5 kpc, Bq - 2 pG and 
Be - 2 pG. and 



Diir) = 



+ 1 r > 1.5 kpc 
-16 kpc < r < 7.5 kpc 
+ 15 kpc < r < 6 kpc 
-1 r < 5 kpc 



(B.5) 



While the halo field is given by 
1 r 



B"(r,z)^B^- 



1 + 



,w\ 



R 



■exp - 



R, 



H\ 



^0 



(B.6) 



and the parameters are Zq = 1.5 kpc, B^ =10 pG, Rq = 
4 kpc, z" = 0.2 kpc (for \z\ < z"), and zf = 0.4 kpc (other- 
wise). 



Appendix B: Magnetic field models 

For the convenience of the reader we present the parameteriza- 
tion s of the galacti c magnetic field folowing Page et aTl (|2007|) 
and lSun et al.l (1200 8). Everything is in the usual cylindrical co- 
ordinates, whose origin lies at the Galactic center. The x-axis of 
the coordinate system points in the opposite direction of the Sun, 
while the z-axis points towards the Galactic north. 

* The final RM and UDM value is defined as the sum of the highest 
possible resolution sub-beam to all parent beam contributions. Thus one 
intensity pixel might have several corresponding RM and UDM pixels, 
since those are computed for higher resolutions, while the former aver- 
ages over the same higher resolution computations. 



